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Abstract 



The site percolation thresholds p c are determined to high precision for 
eight Archimedean (uniform) lattices, by the hull- walk gradient-percolation 
simulation technique, with the results p c = 0.697043, honeycomb or (6 3 ), 
0.807904 (3, 12 2 ), 0.747806 (4,6,12), 0.729 724 (4,8 2 ), 0.579498 (3 4 ,6), 
0.621819 (3,4,6,4), 0.550 213 (3 3 ,4 2 ), and 0.550 806 (3 2 ,4,3,4), and errors 
of about ±2 x 10 -6 . (The remaining Archimedean lattices are the square 
(4 4 ), triangular (3 6 ) and Kagome (3,6,3,6), for which p c is already known 
exactly or to a high degree of accuracy.) The numerical result for the (3, 12 2 ) 
lattice is consistent with the exact value [1 — 2 sin(7r/18)] 1 / 2 , which we also 
derive. The values of p c for all Archimedean lattices are found to be lin- 
early related to the density of sites within an error of about 1%, which is 
more accurate than correlations based solely upon the coordination number. 

Comparison with nonuniform lattices is also made. 
PACS numbers(s): 64.60Ak, 05.70.Jk 
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I. INTRODUCTION 



The Archimedean (Arch) tilings were first elucidated and published by Kepler (a trans- 
lation of which can be found in [l|]). These lattices received their name from references in 
Kepler's paper to Archimedes' descriptions of regular solid polyhedra, which are related to 
these 2d lattices. The significance of the Arch tilings or lattices, which are shown in Fig. [I], 
is that they are the complete set of lattices having infinite tessellation in which all vertices 
are equivalent. (They are also called "uniform" lattices.) This property has made them 
useful in the study of mathematics 0, crystallization P,||], statistical mechanics H, as well 
as percolation ||. The familiar square, triangular, honeycomb and Kagome lattices are all 
Arch. 

The primary goal of this study was to determine precise values of the site percolation 
threshold p c for all Arch lattices for which accurate values were not previously known. A few 
years ago, d'Iribarne, et al. determined p c for all but one of the Arch lattices (which they call 
"mosaics") to about three significant figures ||, using a method based upon the analysis 
of the minimum spanning tree of clusters J7|. More recently van der Marck (vdM) |||| 
determined the thresholds for three of these lattices to nearly four figures. Furthermore, the 
threshold for the honeycomb lattice, a common one for percolation studies, was previously 
known only to about four figures J9JI2",IB| . Here, we extend the precision of these thresholds 
to six significant figures. We did not consider the square, triangular, and Kagome lattices, as 
p c (site) is either known exactly (triangular and Kagome @]), or has already been measured 
to a high degree of precision (square for these cases. 

The second goal of this study was to explore the dependence of p c upon lattice character- 
istics and determine to what extent a universal formula can be found for 2d site percolation. 
Recently there has been renewed interest in this question with the extensive work of Galam 
and Mauger (GM) PfjSJ and vdM @,|J2|^4[. The Arch lattices provide an excellent 
resource for studying this question because they possess a range of coordination numbers 
from 3 to 6 and a variety of other lattice characteristics. 
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A. Nomenclature 



Unfortunately, neither Kepler nor Archimedes (nor anyone else to our knowledge) as- 
signed common names to all the Arch lattices. However, Griinbaum and Shephard have 
introduced a general notation, which we use here, to catagorize such lattices in terms of the 
set of polygons which surround each vertex, (n" 1 , n^ 2 , . . .). Going clockwise around a vertex, 
the numbers rij denote the number of sides of each polygon, and the superscript dj refers to 
the number of these polygons adjacent to each other. For example, the triangular lattice has 
six triangles around a given point, and is designated (3 6 ). The labels for all the Arch lattices 
are shown in Fig. 1. Note that the (3, 12 2 ) lattice has been called the "3-12" lattice in the 
literature, and in places the (4, 8 2 ) has been given the name of "bathroom tile" lattice. 

For convenience in our work we did assign nicknames to the lattices, which are listed in 
Table |. The (4, 6, 12) is called the "cross" lattice because of the cross in Fig. § the (4, 8 2 ) 
lattice is named after our local shopping mall, Briarwood, which (like many such malls) has 
tiles of this shape on its floor; (3 3 ,4 2 ) is called the "direct" lattice because of its directed 
structure, etc. The origins of the rest of the names are, we hope, somewhat more obvious. 



II. METHOD 



A. Hull Gradient Walk 



To find p c we employ the hull-gradient method Pq] , which we previously used to de- 
termine p c for site percolation on the square lattice [2~7| and for bond percolation on the 



Kaeome lattice H|, both to more than six significant digits of precision. 



In gradient percolation |]29| , |30| , a linear gradient of occupied sites is applied in the vertical 



direction of the lattice. As the height increases, the concentration of occupied sites also 
increases. This gradient forms two "land masses" within the lattice: a continuous region 
of occupied sites located at the top of the lattice and a continuous region of vacant sites 
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located at the bottom (Fig. Q). The boundary between these two regions forms a path whose 
average height provides an estimate for p c . 

To make this method efficient for finding p c , we employ a hull-generating walk |3l| . |32|j 
which simultaneously generates and identifies the interfacial boundary. The status (whether 
occupied or vacant) of a site encountered by the walk is determined by generating a random 
number and comparing that number with the occupation probability for that height. If the 
walk does not arrive at a given site, the status of that site will not be determined, and no 
random number will be generated for it, contributing to the efficiency In fact, we believe 
that the hull-gradient method is the most efficient way to determine p c for 2d lattices. It is 
also quite simple to program, as it involves no lists or cluster labelling algorithms. 

For each gradient |Vp|, an estimate of p c may also be obtained by taking the ratio of 
occupied to total (occupied plus vacant perimeter) sites belonging to the hull [pO] . As the 
size of the system is increased or the gradient is decreased, p c (|Vp|) approaches p c linearly 
in |Vp|. Thus a simple extrapolation of the data for finite gradients gives the infinite-system 
value. In practice, we considered gradients sufficiently small (typically of the order of 
Ap/lattice spacing) that the extrapolation from the final point to zero gradient was of the 
order of the error bars of that last point. Additional details of the method are given in 
IIPl. 



B. Computer Techniques 

In the simulation, the status of all points visited by the walk must be stored in computer 
memory. To accomplish this efficiently, all of the Arch lattices were transformed to align on 
rectilinear grids, as shown in Fig. [| Note that the gradient was not applied directly to these 
squared-off forms, as some lattices had to be distorted in the vertical direction to get them 
in this form. Instead, all numerical calculations were performed so that the height of the 
point in the original, undistorted, lattice was used to find p(z). Previously we showed, for 
bond percolation on the Kagome lattice, that assigning p(z) this way yields the best linear 
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scaling of p c vs. gradient [28[|. For two lattices, we considered two possible orientations of 
the lattice and found similar results, as discussed in Appendix A. 

The lattices were initialized by filling the top half of the first column with occupied sites 
and the lower half with vacant sites, which prevents the walk from closing on itself at the 
start of the algorithm (Fig. |2[). Periodic boundary conditions were applied in the horizontal 
direction, and each new column to the right was cleared as it was visited. This allowed the 
simulation to run indefinitely and have essentially no boundary effects from the horizontal 
ends of the system. The maximum distance the walk traveled horizontally from the front 
was tracked to detect if wraparound errors occurred; if they did, the system size and/or 
gradient was adjusted accordingly, and the run was restarted with the expanded system. 

Random numbers were generated using the shift-register sequence generator R7(9689) 
[ 3|,|3l defined by: 



Xn — Xn-471 A ^n-1586 A ^n-6988 A ^n-9689 j (1) 

where A is the bitwise exclusive-or operation. This 'four tap' generator is equivalent to 
decimating by 7 (taking every seventh term) of the sequence generated by the two-tap rule 



R(9689) (x n = x n _47i A x n _9 68 g) given by Zierler f3"5|l , where the decimation has the effect 
of greatly reducing the three- and four-point correlations of the two-tap generator ||33|| . 
Our previous work has shown that this generator does not appear to introduce errors in 
simulations of this kind IE 



III. RESULTS AND DISCUSSION 



A. Percolation thresholds 



For each lattice, p c was plotted as a function of the gradient. Fig. f| shows a representative 
plot for the (3, 12 2 ) lattice. In all linear relationship was observed between the 

magnitude of the gradient and the estimate p c (|Vp|), as found previously for the square 
and Kagome lattices p8| . The value at the intercept of the |/-axis represents p c for an 



infinite lattice (zero gradient), and these values are reported in Table |. Also shown in Table 
| is the quantity of random numbers generated for each simulation, which is identical to 
the total number of occupied and vacant sites belonging to the hull, because one random 
number is generated for each new site visited. These simulations consumed several months 
of computer time on Sun and HP workstations. 

For all our values, the error (which is purely statistical in origin) is of order 10~ 6 . Our 
results are consistent with d'Iribarne et al.'s three-digit values, if we assume an error of 
±0.003 in their data. Our result 0.697043(2) for the honeycomb lattice is consistent with 
the previous values 0.6973(8) of Vicsek and Kertesz fl~2]j , 0.6962(6) of Djordjevic et al. fll3"f , 
and 0.6971(2) of vdM 0. Likewise, our results agree with the determinations 0.5504(2) for 
(3 3 ,4 2 ) and 0.7298(1) for (4, 8 2 ) made recently by vdM . 



In another recent work, Prea |36j has calculated the distance sequences c n , defined as the 



number of sites n steps (chemical distance) from the origin, for the Arch lattices. Prea found 
a monotonic relation between p c (site) and c = lim^oo inf (c n /n), although three different 
lattices with different p c — the square, Kagome, and (3, 4, 6, 4) — share the same value 
c = 4. This apparent monotonic behavior led Prea to conjecture that p c for the (3 4 , 6) 
lattice (the one Arch lattice not studied previously) is within in interval (0.55,0.6). Our 
result 0.579 498(2) confirms this conjecture. 

B. Exact percolation threshold for the (3, 12 2 ) lattice 

The threshold for the (3, 12 2 ) can be derived exactly. The proof is based upon the 
similarity of structure between this lattice and the Kagome lattice. In Fig. |5], the (3, 12 2 ) 
lattice is shown with the triangle bonds in bold and the bonds between the triangles as 
dashed lines. If the dashed lines are reduced in length, the (3, 12 2 ) lattice transforms into 
the Kagome lattice. Similarly, the sites on both ends of the dashed line in the (3, 12 2 ) 
(marked by the large circle in the figure), taken together, are equivalent to a single point on 
the Kagome lattice. 
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Let both lattices be at the percolation threshold and p be the probability that a given 
site on the (3, 12 2 ) lattice is occupied. Then, p 2 is the probability that both sites at the 
ends of the dashed bonds are occupied, and the bond will permit flow. The dashed bond is 
equivalent to a single point on the Kagome. At the percolation threshold, the probability 
that a site is occupied on the Kagome lattice is 1 — 2sin(7r/18). Therefore, for the (3, 12 2 ) 
lattice, 

p c = [1 -2 sin(7r/18)] 1/2 = 0.807900 764... . (2) 

Our numerical estimate 0.807904 is consistent with this result within the error bars 
±0.000 003. This agreement provides a further validation of both our algorithm and er- 
ror analysis. 

C. Galam and Mauger's universal formula for p c 

Early in the work of percolation, thresholds were related with coordination number q 
using relatively simple formulas. More recently, GM have attempted to refine this by 

introducing different classes of lattices, characterized by the constants po, a and b in their 
formula 

p c =p [(d-l)(q-l)]- a d\ (3) 

where d is the dimensionality. By defining two different classes, each with distinct sets 
of constants po and a for site and bond percolation (with b = a for bond and b = for 
site percolation), GM were able to obtain a fairly good fit of p c for a number of lattices of 
various dimensionality, with a fifth set of constants po and a required for systems of very 
high dimensionality. Their class 1 contains all 2d lattices except the Kagome, while class 2 
comprises all higher-dimensional systems plus the 2d Kagome. 

Subsequently, vdM [||,^2j considered a number of additional lattices and argued that a 
formula of this type, depending only upon d and q, cannot be very accurate because of the 
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existence of different lattices with differing p c but identical d and q. GM |19j responded that 
their formula is applicable only for "isotropic" lattices in which each site has the identical 
q. (In this present paper, we describe these lattices as "uniform," by which we mean that 
each vertex is identical in terms of its surrounding polygons, allowing for rotations and 
reflections.) For other lattices, GM reformulated @ by reinterpreting q to be an effective 
coordination number q e g (not necessarily the average value q). Because q e s cannot be found 
independently, this reinterpretation has the effect of turning their formula into a correlation 
between site and bond thresholds for a given lattice rather than a more general correlation 
for p c , which is developed further in ]2l| and discussed in PJII]. However, because we do 
not consider thresholds for bond percolation here, we do not address this aspect of GM's 
work in this paper. 

In Fig. ||we plot our results using the same axes as GM [p0[ , ln(l/p c ) vs. ln(g— 1). We also 
exhibit GM's two fitting formulas for site percolation, which are straight lines on this plot. 
Clearly, many of the Arch lattices do not fall near either of these two lines. Evidently, GM's 
approach (|3|) cannot accommodate these additional uniform lattices, unless one introduces 
still more lattice classes, which is counter to the concept of universality. 



D. A new fitting formula for p c 

For the Arch lattices, one can observe the trend that p c increases with increasing "open- 
ness" of the lattices. To characterize the latter, we introduce the quantity p equal to the 
density or number of lattice sites per unit area, assuming that all bond lengths are unity. 
Its inverse, the area per site, can be determined by drawing lines which bisect the centers of 
each polygon surrounding a given site, as in a dual-lattice construction, and summing the 
enclosed area. The tiles defined by these bisector lines clearly fill the entire lattice. The 
total area per site is simply 1/3 the area of each triangle, 1/4 the area of each square, etc., 
surrounding that site. For Arch lattices with uniform bond lengths, all polygons are regular. 
The area of 1/n-th of a regular n-gon, with unit edge length, is given by 
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. 1 7T 

A n = - COt - , 

4 n 



(4) 



with some numerical values given in Table [TT]. For a lattice characterized by vertices 
(n® 1 , n^ 2 , . . .), we then have 



n -l 



£Jj cot — 



(5) 



The resulting values of p for the Arch lattices are listed in Table |TJ. In Fig. |7| we show 
a plot of p c vs. p, and indeed one can see that the correlation between p c and p is good. 
Furthermore, that correlation can be well represented by a simple linear function. A least- 
squares fit yields 



Pc 



1.0405 - 0.4573p . 



(6) 



The errors involved in using this formula are listed in the column A in Table PTT| , and are 
generally within about 0.01, in contrast to errors up to 0.05 if GM's correlations were used 
(taking the best class in each case). In Fig. |8] we plot the errors from using these different 
fitting formulas for the various lattices, clearly showing how GM's formulas are each accurate 
for certain lattices only, while (||) works fairly well for all of the lattices. The rms errors 
from using @ for the 11 Arch lattices equals 0.0075. 

It turns out that our fit (f|||) is related to a correlation for p c given by Scher and Zallen 
nearly three decades ago [p9[| . These authors introduced a "filling factor" / defined as the 
fraction of space occupied by disks of radius 1/2 placed at each lattice site, and found for 
four lattices in 2d, 



fp c « const. = 0.44 ± 0.02 . 



(7) 



Now, / is simply 7rp/4, so that their correlation also relates to the density of sites. However, 
while their hyperbolic relation (|7|) provided a good fit to lattices they considered (to the 
accuracy available at that time), it does not capture the behavior for the other lattices, 
especially those of low p, as shown in Fig. ^, where (|7|) is also plotted. 



Fairly recently, d'Iribarne et al. |J considered a correlation for p c of the Arch lattices 
that is also closely related to p or /. They compared p c to a variable m which represents 
the average edge length of a minimally spanning tree on a complete lattice, normalized by 
the area. Because on the complete graph there is one edge per site, m is identical to p 1 / 2 . 
Their plot of p c vs. m correlates the data to a single curve just as our Fig. [7| does, however, 
in terms of m the behavior is quite curved, which d'Iribarne et al. fit to the quadratic 

p c = 0.685 + 0.799m - 0.899m 2 . (8) 

In fact, this fit reduces the rms error of our results over the 11 lattices to 0.0057. (Note, only 
a slight improvement is obtained by adjusting these coefficients for our new values of p c .) 
However, the expense of this improvement is an extra parameter in (|8]) over (|J). Indeed, 
one can also extend fl6|) to a quadratic in p, with the result 

p c = 0.9472 - 0.2181p - 0.1439p 2 , (9) 

which gives a comparable rms error over the 11 lattices, 0.0054, but again at the expense of 
a more complicated formula than (P). 

Along these lines, it is interesting to note that an excellent fit can be achieved using 
results from only the three lattices where p c is known exactly (which happen to span the 
whole range of p c values). That fit, 

p c = 0.9466 - 0.1972p - 0.1642p 2 , (10) 

yields an rms error over the 11 Arch lattices of 0.0062, only slightly worse than (|Sp or ([5|). 
This fitting formula depends upon no adjustable parameters nor any Monte- Carlo- measured 
values. 

E. Extension to non-uniform lattices 

We have found that (|§-^) may also be applied to non-uniform lattices, if in calculating p 
by ([D, a n is interpreted as the average number of polygons of type n over all vertices. For 
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the non-uniform lattices (also referred to as anisotropic |2(J), many of the polygons are not 
regular, and it is not always possible to make all bonds of the same length, which were the 
two assumptions in deriving (|5]). Thus, we now consider ([5]) to represent a topological or 
graphical weight rather than an actual geometric property of the lattices. 

We looked at seven non- uniform lattices considered by vdM [^,[| and/or GM [20| and 



whose values of p c are known (we did not determine any of these p c 's here): the bowtie, 
bowtie dual, pentagonal, (4, 8 2 )-matching, dice, Penrose, and Penrose dual lattices. The 
first five of these lattices are shown in Fig. |9|, while the quasicrystalline Penrose and Penrose 
dual are shown for example in |38| and |H| . The form of the Penrose lattice used here is 
the rhomb or P3 form |[25|| . 

For the dice lattice, | of the vertices are (4 6 ) while | of them are (4 3 ), so on the average 
the vertices are (4 4 ). Likewise for the (4, 8 2 )-matching, we have |(3 2 ,8 2 ) + |(3 2 ,4,8) = 
(3 2 ,4 2 / 3 ,8 4 / 3 ), for the bowtie |(3 4 , 4 2 ) + ±(3 2 ,4 2 ) = (3 3 ,4 2 ), for the bowtie dual |(4 2 ,6 2 ) + 
|(4, 6 2 ) = (4 4 / 3 ,6 2 ), and for the pentagonal |(5 4 ) + |(5 3 ) = (5 10 / 3 ). For the Penrose 
lattice, the average vertices are simply (4 4 ), while for the Penrose dual lattice the average 
vertices are Q (3 a3 , 4 a4 , 5 as , 6 ae , 7 ar ) with a 3 = 3(x 2 + x 4 ) = 3(3 - Ax) = 1.583592135, 
a 4 = 4x 5 = 4(5x-3) = 0.360679775, a 5 = 5(x 3 + x 6 ) = 5(4- 6x) = 1.458980338, a 6 = 6x 7 = 
6(13x - 8) = 0.206651122, and a 7 = 7x 6 = 7(5 - 8x) = 0.39009663, where x = (y/E-l)/2 
is the inverse of the golden ratio. Note that in general a n satisfy Xm(l — ^-/ n ) a n — 1 and 
J2 n a n = Q, where q is the average coordination number. 



These average vertex numbers are listed in Table fV|, along with p c values, p values, p^ st ' 
and error A from using (|5|-|6|). Interestingly, the A for all 7 of these lattices are generally 
smaller that those for the Arch lattices — with an rms value of 0.0046 vs. 0.0075. (Note 
that some of these A's are however close to the precision of the value of p c .) Surprisingly, 
even the Penrose dual lattice, which incorporates polygons with from 3 to 7 sides and a very 
irregular structure, is also well modelled by this formula. Evidently, this approach, which is 
based upon p given by (^), appears to be quite robust. 

One could conceivably use these additional data points from the non-uniform lattices to 
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refine the fit of the formulas for p c {p) in Section |I [D|. However, the improvement turns 



out to be marginal, and considering the relatively low precision of the values of p c for the 
non-uniform lattices, the justification for doing this is doubtful. 

Note that if one retained the definition of p literally as the density of sites for these 
lattices, rather than using the definition above, then the fit of the non-uniform lattices with 
(g) would be very poor. An example is given by the dice lattice, which has a site density 
identical to that of the triangular lattice, 2y / 3/3 = 1.1547. The p c for this lattice, 0.5851 
||, is substantially above the value for the triangular, 1/2, and the data point for this case 
would fall well above the line in Fig. [f[ 

It is also interesting to note that we have now considered three different lattices with the 
same average vertex environment of three triangles and two squares: the two Arch lattices 
(3 3 ,4 2 ) and (3 2 ,4, 3,4), and the non-uniform bowtie lattice, which have p c 's of 0.550 806, 
0.550 213, and 0.5475(8) respectively. The closeness of these values supports the idea that 
p c is determined principally by the local nature of the lattice as characterized by the vertex 
numbers. For the two Arch lattices, the p c 's are nearly identical, and close to the prediction 
by (|6|) of 0.5503. The bowtie lattice, which has the same vertex numbers as the others on 
the average only, has a value of p c that is slightly lower, by about 0.003. 

The inadequacy of using q alone to estimate p c is also apparent, as the (3 4 , 6) lattices 
has the same q = 5 as these three above, but a p c that is substantially higher, 0.579 498. 
The results for the non-uniform lattices are also included in Fig. using q for q. 

A similar situation occurs for the square (p c = 0.592 746), dice (p c = 0.5848), and Penrose 
(p c = 0.5837) lattices, whose (average) vertices are (4 4 ) in all cases. Again, the non-uniform 
lattices have a slightly lower p c than the uniform one with the same vertices, with the "most" 
non-uniform one (the Penrose) having the lowest p c . And again, other lattices with q = 4 
but different average vertex numbers (there are four of them here) all have much different 
(and higher) p c 's. 

van der Marck has pointed out p3] , pi][ ] that adding an extra site at the center of any 



triangle on a 2d lattice does not change the value of p c (site). If these extra sites are 
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included in the calculation of p, then p would change and the predictions of |6] would be 
poor. Therefore, in applying these formulas, one must disregard these superfluous sites. 
However, there are undoubtedly other classes of lattices where |5|-|6] will not work well — and 
will not be so easily fixed. 

F. Further variations on the fitting formula 

Although a precise universal formula for p c based upon the vertex numbers alone is 
impossible, it is interesting to explore variations to (HH) to see what improvements can be 
made. According to (|5|), the variable p~ l is an average over the nearest neighbors polygons 
weighted by A n . To see how close this weighting compares to an optimal one, we carried out 
a regression analysis, allowing all A n as well as the linear fit of p c {p) to vary, and minimizing 
the error from the linear fit using data from all 18 lattices given above, both Arch and 
non-uniform. This procedure yielded the A' n given in Table [II] and the fit 

p c = 1.1383 - 0.5503p' , (11) 

where p' = J2 a nA' n . Note that we kept A' A = 0.25, as one weight can be fixed arbitrarily. 

The errors using this fit are generally lower than those that result from (^); in fact the 
rms error over all 18 lattices decreases substantially from 0.0065 when @ is used to 0.0028 
when (|ll|) is used. However, this formula is clearly less predictive than eq. which had 
no adjustable parameters in the definition of p. 

The A' n listed in Table [II] fall close to, but in a somewhat narrower range than, the A n . 
(Note A' 7 is somewhat out of line with the rest, being bigger than A?, but was derived from 
only one lattice, the Penrose dual.) These results suggest that there might be a topological 
weighting rule that is better than the effective area per site, ([5]). However, we do not know 
what that weighting rule may be. 

To test if a non-linear relation between p c and p' may be better, we carried out the same 
regression analysis as above, but allowed for a quadratic p c (p)- Interestingly, the A' n came 
out nearly the same, and the behavior of p c {p) was very nearly linear. 
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Using d'Iribarne et al.'s form with m = y/p* as the variable and assuming a quadratic 
relation between p c and m, we found optimal weights A' n now much closer to the A n . These 
various and partly contradictory results suggest that, while some improvements can be made, 
our original assumptions — that A n is the actual polygon area, and that p c {p) is linear - 
are not unreasonable. 

Of course, there are deeper assumptions to this whole approach — that p~ l is a linear 
function of the a n , and that p c is a simple linear or quadratic function of p — which we have 
not explored. And ultimately, to get a more precise formula for p c , one must also take into 
account the actual arrangement of the polygons around each vertex, not just the number of 
each type of polygon. Clearly, further work can be done in understanding the relation of p c 
to lattice type. 



G. Conclusions 

We have determined the p c (site) for the honeycomb and seven less-common Arch lattices 
to nearly six significant figures, about 100 to 1000 times more precise than previously known 
- except for the (3 4 , 6) lattice, whose p c was never determined previously. The value for 



the latter lattice is within the conjectured bounds of Prea ||36|| . We find a fairly good linear 



correlation between p c and p (and related to the correlations of Scher and Zallen [39] and 
of d'Iribarne et al. 0), although its accuracy is far below what can be found in simulation. 
This correlation can be extended to non-uniform lattices, and can be improved somewhat 
if the weightings A n are adjusted to the A' n . The correlation with p is substantially more 
accurate than correlations based primarily upon coordination number, such as those of GM 



All our results here are restricted to 2d and to site percolation. It would be interesting 
to see if these ideas on relating p c to the local graphical environment can be extended to 
bond percolation and to higher dimensionality. It would also be interesting to see how well 
GM's ideas relating p c (site) and p c (bond) on the same lattices hold for the Arch lattices. 
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APPENDIX A: IMPORTANCE OF GRADIENT ORIENTATION 

In order to verify that the direction that the gradient is applied on a lattice has no 
effect on p c , we calculated p c for both the honeycomb and (3 3 ,4 2 ) lattices after they had 



been rotated 90 degrees. Fig. [10] shows the original orientation of the simulation and the 
orientation after the lattice has been rotated 90 degrees. The p c for the honeycomb lattice 
was measured as 0.697 043±0.000 002 and then 0.697 046±0.000 003 after 90 degree rotation. 
Likewise, for the (3 3 , 4 2 ) lattice, p c was measured as 0.550 213±0.000 002 and then 0.550 21 1± 
0.000 002 after 90 degree rotation. In both cases, the p c values were identical within the error 
of the method, verifying that the orientation of the lattice has no effect on determining p c 
in our method. 

^Electronic mail: rziff@engin.umich.edu 
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FIGURES 

FIG. 1. The 11 Archimedean (Arch) lattices, in which all vertices are equivalent. Lattices are 
designated using the notation of Griinbaum and Shephard [23] as explained in the text. 

FIG. 2. The hull-generating walk along a percolation gradient for the honeycomb lattice. 
Filled circles denote occupied sites and heavy lines show the hull of the percolating region. The 
arrow points to the starting point of the walk. 

FIG. 3. The Arch lattices transformed to a square array, for use in the computer simulation. 
Some lattices were distorted in the vertical direction, but in all cases the actual lattice height, 
rather than the height in the square array, was used to determine the site occupation probability. 

FIG. 4. Percolation threshold vs. lattice gradient for the (3, 12 2 ) lattice. As the gradient 
becomes smaller, the lattice approaches infinite size. The percolation threshold can be estimated 
from the y-intercept of a linear regression of results from finite lattices. 

FIG. 5. Derivation of the exact percolation threshold for the (3, 12 2 ) lattice. Replacing the 
dashed bonds with single points leads to the Kagome lattice. 

FIG. 6. Thresholds plotted as In l/p c vs In {q — 1). (•) regular lattices (from left to right: 
honeycomb, square and triangular), (□) Kagome, (A) less-common Archimedean lattices, (x) 
non-uniform lattices. Dashed line: GM's class 1 formula; solid line, GM's class 2 formula. Evi- 
dently, many lattices do not conform to GM's classification scheme. 

FIG. 7. Percolation thresholds as a function of p, with the same symbols as in Fig. |6|. The 
linear fit is given in eq. (@). Also shown are is correlation of Scher and Zallen, eq. (0) (dashed line). 



FIG. 8. Error from using GM's class 1 formula (□), class 2 formula (A), and our fit (||-^) (•), 
for the 18 lattices, plotted sequentially in the order they are listed in Tables [LTl| and [lV| respectively 
Formulas (^), (||) or (11) improves the fit of (|§-§) even further. 
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FIG. 9. Five of the non- uniform lattices used to test the fitting formulas. 

FIG. 10. Two orientations for the (3 3 ,4 2 ) and honeycomb lattices used to test the influence of 
lattice orientation in our method of determining p c . 
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TABLES 




TABLE I. Measured values 


of p c for site percolation 


on the Archimedian lattices. Lattice 


nicknames in quotes are ours. 






Lattice 


Measured p c 


Total Random Numbers Generated 


(3, 12 2 ) "star" 


0.807904(3) 


4.1 x 10 11 


(4,6,12) "cross" 


0.747806(3) 


2.6 x 10 11 


(4,8 2 ) "Briarwood" 


0.729 724(2) 


2.6 x 10 11 


(6 3 ) honeycomb 


0.697043(2) 


4.0 x 10 11 


(3,4,6,4) "bounce" 


0.621819(2) 


2.9 x 10 11 


(3 4 ,6) "bridge" 


0.579498(2) 


2.9 x 10 11 


(3 2 ,4,3,4) "puzzle" 


0.550 806(2) 


2.8 x 10 11 


(3 3 ,4 2 ) "direct" 


0.550 213(2) 


2.9 x 10 11 


TABLE II. Weights A n @, 


and optimized weights A' n 


which yield the best linear fit of p c (p') 


over all 18 lattices. 






n 


A n 


A' 


3 


0.144338 


0.14400 


4 


0.25 


0.25 


5 


0.344095 


0.33640 


6 


0.433013 


0.41653 


7 


0.519130 


0.52539 


8 


0.603553 


0.55041 


12 


0.933013 


0.75490 
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TABLE III. Values of p c and lattice characteristics for all 11 Arch lattices. A is the error when 
@ is used to fit the data. 



Lattice 


Pc 


Ref. 


q 


pby © 


pf" by (|) 


A 


(3,12 2 ) 


0.807904 




3 


0.4974 


0.8130 


-0.005 


(4,6,12) 


0.747806 




3 


0.6188 


0.7575 


-0.010 


(4,8 2 ) 


0.729 724 




3 


0.6863 


0.7266 


0.003 


(6 3 ) honeycomb 


0.697043 




3 


0.7698 


0.6885 


0.009 


(3, 6, 3, 6) Kagome 


0.652 703 6 





4 


0.8660 


0.6444 


0.008 


(3,4,6,4) 


0.621819 




4 


0.9282 


0.6160 


0.006 


(4 ) square 


0.592 746 




4 


1 


0.5832 


0.010 


(3 4 ,6) 


0.579498 




5 


0.9897 


0.5879 


-0.008 


(3 2 ,4,3,4) 


0.550 806 




5 


1.0718 


0.5503 


0.0005 


(3 3 ,4 2 ) 


0.550 213 




5 


1.0718 


0.5503 


-0.0001 


(3 6 ) triangular 


0.5 


[14] 


6 


1.1546 


0.5124 


-0.012 



TABLE IV. Data for non-uniform lattices. Here vertex numbers represent averages. 



Lattice 



Ref. 



P by 



A 



(4,8 2 )-matching (3 2 ,4§,8l) 0.6768(2) 

bowtie dual (43, 6 2 ) 0.6649(2) 

pentagonal (5T) 0.6476(2) 

Penrose dual (see text) 0.6381(3) 

dice (4 4 ) 0.5848(2) 

Penrose (4 4 ) 0.5837(3) 

bowtie ex (3 3 ,4 2 ) 0.5474(2) 



IE3 
III 
ill 



§0 



4 
4 
4 
5 



0.7936 
0.8338 
0.8719 
0.8987 

1 

1 

1.0718 



0.6776 
0.6592 
0.6418 
0.6295 
0.5832 
0.5832 
0.5503 



-0.0008 
0.006 
0.006 
0.009 
0.002 
0.0005 
-0.003 
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